Measuring our Peculiar Velocity by "Pre-deboosting" the CMB 
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It was recently shown that our peculiar velocity /3 with respect to the CMB induces mixing 
among multipoles and off-diagonal correlations at all scales which can be used as a measurement 
of f3, which is independent of the standard measurement using the CMB temperature dipole. The 
proposed techniques rely however on a perturbative expansion which breaks down for i > 1//3 « 800. 
Here we propose a technique which consists of deboosting the CMB temperature in the time-ordered 
data and show that it extends the validity of the perturbation analysis multipoles up to I ~ 10000. 
We also obtain accurate fitting functions for the mixing between multipoles valid in a full non- 
linear treatment. Finally we forecast the achievable precision with which these correlations can be 
measured in a number of current and future CMB missions. We show that Planck could measure 
the velocity with a precision of around 60 km/s, ACTPol in 4 years around 40 km/s, while proposed 
future experiments could further shrink this error bar by over a factor of around 2. 



I. INTRODUCTION 

The dipole of the CMB is measured to be much larger 
than the other multipoles and this is usually attributed 
to a Doppler effect due to our peculiar velocity /3 with 
respect to the CMB rest frame. Under this assumption 
we can infer, by combining the measured temperatures 
of WMAP dipole [1] with the COBE monopole [2, 3], its 
direction to be I = 263.99° ± 0.14°, b = 48.26° ± 0.03° 
in galactic coordinates, and its modulus to be (3 = \f3\ = 
(1.231 ± 0.003) x 10~ 3 . These very precise numbers rely 
however entirely on the above assumption, but generi- 
cally the CMB dipole is not necessarily due only to a 
relative velocity. One way to test this assumption was 
proposed in [4] where it was shown that our peculiar ve- 
locity could also be measured using the asymmetry in 
the location of the peaks of the power spectrum between 
forward and backward hemispheres, achieving a possible 
detection of /3 at 2 — 3cr for resolution £ = 1000 — 1500. 
Recently it has also been pointed out [5-7] that all the 
CMB multipoles at m have a correction due to our local 
peculiar velocity because the primordial anisotropies are 
distorted by the Doppler and aberration effects. This 
shows up as a correlation between different multipoles I 
and could be used as an alternative way of measuring our 
velocity, as an independent consistency check. This fact 
might offer also a way to test the isotropy of the Uni- 
verse on very large scales: it is known in fact that the 
CMB sky and other observations seem to exhibit some 
anomalies on the very large scales and this effect offers 
an observational handle which could either confirm and 
make more robust the standard assumptions or perhaps 
point to global anisotropies of the Universe. 

The analysis of these CMB correlations has been per- 
formed [5-8] relying upon a Taylor expansion in orders of 
P of boost effects (Doppler and aberration) on the multi- 
pole coefficients af m , where X stands for T (temperature) 
or E, B (the E and B modes of polarization). It turns 



out, however, that each of the series coefficients brings 
together ever higher powers of £, effectively transform- 
ing the expansion into one of powers of the product (3 1. 
Since /3 = 1.231 x 10~ 3 this means that for i > 800 the 
series can no longer be relied upon. We shall therefore 
distinguish between two regimes of angular scales where 
the the above correlations show up. In the first regime on 
large scales, I -C 1//? — 800, we can use a perturbative 
approach and therefore precisely predict that there are 
correlations of (D((3£) between neighbor ts and look for a 
signal in the CMB in the form of correlators a* im at+i m i . 
Moreover it was shown in [7] that one can non-trivially 
and very efficiently define three different estimators (for 
m! — m = 0, ±1) in order to measure directly the three 
cartesian components of (3, without having to scan all 
the possible directions in the sky (i.e., without the need 
to compute and minimize a numerical \ 2 for the correla- 
tions for each direction in a grid of {9, cj>} coordinates). 

The multipoles in the complementary regime (£ > 1//3) 
cannot be so easily treated but nevertheless also carry in- 
formation about our peculiar velocity in the form of sim- 
ilar correlations, and it turns out that such information 
is needed in order to reach a signal-to-noise ratio larger 
than 1. In this regime the deviation angle 9 — 9' due 
to aberration is larger than the angular scale of inter- 
est l/£, so things are much more complicated. There 
are in fact correlations which are large (0(1)!) and 
nonzero also for distant ^'s given by a transformation 

<4m = X^=o K# t m <4^L > between the aberrated [A] 
frame and the primordial [P] correlations. One in princi- 
ple would have to compute the matrix elements K^i m , 
sometimes referred to as the aberration kernel, which are 
integrals of spherical harmonics with different arguments, 
plug the matrix elements into all the possible two-point 
correlation functions and compare with the data. In fact, 
as we will show in Section II the correlations in non- 
neighboring multipoles (a* lm ai± nm i , n > 1) also carry a 
measurable signal, which should be taken into account to 
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measure the velocity with reasonable precision. While in 
principle straightforward this procedure has some disad- 
vantages: (i) this can be a heavy and delicate numerical 
task, because of the highly oscillating integrands and also 
because of the huge number of correlators that one would 
have to consider (future experiments propose to measure 
all multipoles as far as £ ~ 3000+ for both temperature 
and polarization); (ii) it is not obvious to understand 
in this case whether the three simple estimators can be 
written explicitly for the three cartesian components of 
the velocity for any n in a^ m ag± n m > , so perhaps the pro- 
cedure would have to be carried out scanning the sky in 
all possible directions which would probably make this 
approach even more expensive in terms of computational 
time. 

Even though we explore this approach further and pro- 
pose a solution (in Section II) to the first of the two chal- 
lenges listed above, the main point of this paper is to sug- 
gest a very simple trick which may be used to overcome 
all of these problems and measure directly the three com- 
ponents of (3 from a map in a much faster and straight- 
forward way. Using this result we are able to predict 
with which precision we can measure (3 for several future 
experiments. This trick relies on the fact that we may 
already assume to know the direction and modulus of the 
velocity (3 from the CMB dipole and that we can the use 
the other multipoles with their correlations as an inde- 
pendent consistency check, in order to confirm (or not) 
the assumed value of (3 up to some precision. 

The trick (described in more detail in Section III) is 
as follows. Given the central assumed value of the ve- 
locity /3 dip = /3 d * p ± <5/3 dip indicated by the dipole, we 
may take the CMB map T(8, <fi) and "deboost" it by a 
Lorentz transformation into the frame with opposite ve- 
locity — /3 d j p . As discussed in [9] such transformation 
should be carried in the CMB time-ordered data (TOD), 
i.e., before data treatment to extract the harmonic mul- 
tipole coefficients ai m and thus before one constructs 
the CMB temperature maps. Thus we call this tech- 
nique "pre-deboosting" the CMB data. In the new (pre- 
deboosted) frame the residual velocity /3 rcs compared 
to the CMB frame is expected to be given just by the 
residual error Sf3 dip on the experimental determination 
of /3 dip . As quoted above, such error is currently ap- 
proximately |#/3 dip | — 3 x 10~ 6 . As a consequence 
the correlations due to Doppler and aberration in this 
new frame are expected to be of the order 0(j3 Tes £) for 
£ <C l/(f3 res ) ~ 3 x 10 5 and so we could safely use the first 
order perturbative equations to compute the correlation 
functions a* lm ag+\ m i up to, say, £ ~ 10000 which is more 
than enough for all future experiments. Note that this 
allows us to use directly again the three efficient define 
estimators (for m! — m = 0,±1) for measuring direction 
and modulus of the residual velocity. 

An unexpected but very interesting prospect is the case 
in which the measured velocity turns out to be different 
from the expected /3 rcs . This would be a clear indication 
that the CMB dipole is not completely (nor around 99%, 



as sometimes stated) due to our peculiar velocity alone. 
Therefore, this would imply that there are other contri- 
butions to the CMB dipole, and it would be interesting 
to understand whether such correlations may distinguish 
even the nature of such contributions: adiabatic, isocur- 
vature perturbations, dipolar lensing or other more exotic 
contributions. 

The ability to measure exotic contributions to the 
dipole is of great interest to test Cosmology on very large 
scales, which could hide non-trivial phenomena, as sug- 
gested by some reported anomalies on the \aw-£ CMB 
multipoles itself [10, 11]. For instance, it could provide 
valuable information about some proposed tilted cosmo- 
logical models in which the dipole arises partly due to pri- 
mordial superhorizon-scale isocurvature fluctuations [12- 
16], which could provide a possible explanation some re- 
cent controversial claims of a high galaxy cluster [17-19] 
and galaxy [20] bulk flow on large scales and it could be 
used as a test of non-standard cosmological models. In 
any case, if /3 res 7^ we would also be able to measure 
the direction of our velocity with an error 86 = 5j3/j3 Tes 
(see [7]). 

This paper is organized as follows. We start by dis- 
cussing in Section II the full non-linear approach to es- 
timate our velocity, and derive a very accurate fitting 
function for the oscillating integrals. We then discuss in 
Section III the pre-deboost technique originally proposed 
by [9] . In Section IV we forecast the sensitivity expected 
from present and future experiments. Finally, we draw 
our conclusions and summarize our results. 



II. THE FULL NON-LINEAR FIT TECHNIQUE 

A. Fitting functions for the aberration kernel 

It was shown in [5-7] that when subjected to an aber- 
ration effect the a^ m coefficients of the spherical harmonic 
decomposition of the temperature (and polarization) con- 
trast transform in the following way 
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where the superscript X stands for either temperature 
(T) or one of the two independent modes of polarization 
[E and B) and where [A] denotes the aberrated coeffi- 
cients, to be contrasted with [P], the primordial (non- 
aberrated) ones. 

In the case of temperature, the exact coefficients of (1) 
are given by [6] (we here follow [7] for the convention of 
the sense of (3 which results in an overall sign change in 
the velocity (3) 
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where 7 = — j3 2 is the standard Lorentz factor 

and 



pm (x) = / 2^+l(^-m)! 



(3) 



and where P™ are the associated Legendre polynomials. 
For polarization the formulae are similar if one makes 
use of spin-weighted spherical harmonics 1 ( S Y™). Fol- 
lowing [5], we get 
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and where in turn s P™(x) = V27T s Y^ m (a;, <^> = 0), the 
spin-weighted spherical harmonics evaluated at <fi = 0, 
which can be written as [231 
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where the sum in r is to be carried out in the range 
max{0, 771 — s} < r < min{£ — s, £ + to}. Note that (5) 
and (6) reduce respectively to (2) and (3) for s = 0. 

The oscillatory nature of these integrals poses a nu- 
merical challenge, which make their direct computation 
very slow for high £. In a recent work [24] a recursive 
method was developed which is claimed to allow fast and 
accurate evaluation of these integrals. In this section we 
follow instead a different route and compute some of the 
integrals just by numerical integration (see Appendix A 
for more details) . Surprisingly, we found out numerically 
that these integrals can be fit very precisely by Bessel 
functions, which greatly simplifies the analysis. 

Following [5] it is convenient to define the quantities 
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where again s represents a spin weight which is for 
temperature and 2 for the E and B modes of polariza- 
tion. Note that we always have < s Gi m < 1/2, the 
lower limit being achieved when \m\ = £ and the higher 



1 For a discussion on spin- weighted spherical harmonics see [21, 22] 



one when to = 0. An exquisite fit for small scales to the 
full non-linear integral for general 777, is given by 
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where J\ is the Bessel function of the first kind. More- 
over, we find that similar relations to the above one ap- 
plies also for non- neighboring correlations (i.e., between 
any £ and £ ± n, n > 1) and also for polarization: 
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which we find to be accurate to around 0.2% for all val- 
ues of £ and to and all values of n. 2 For 77 = one 
cannot apply (9) directly, but we find that an analogous 
fit is given by (with the same precision) 
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As can be seen from the above relations, the tempera- 
ture and polarization aberration kernels become almost 
indistinguishable for small scales (usually for any £ 3> 1, 
or when |m| ~ £, for £ ^> 5). 

A Taylor expansion of the Bessel functions yields to 
leading order 
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which is valid for positive integer values of 77. Expanding 
the above fits in orders of /3 we find that the coefficients 
with the Doppler correction exactly up to leading order 
for any value of n for both temperature and polariza- 
tion. In other words, for a given n the fits are exact to 
order O (/?") . This was confirmed through direct analytic 
integration of the leading order of the Taylor expansion 
in f3 of (2) for specific values of {£, £', m}. This is an in- 
teresting result, as the {£, £ ± 3}, {£, £±4} and so forth 



2 We explicitly checked the above fits for n = {0, 1, 2, 3, 4} and 
£ < 700 for the case of temperature correlations. For polar- 
ization we only tested explicitly n = {0, 1, 2} and i < 100. 
Nevertheless, all results indicate that in both cases the precision 
would remain the same at higher n and at smaller scales. 
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Figure 1. Some of the nonlinear coefficients in (9) for the cases 
n = 1 (solid lines) and n — 2 (dashed) and for different val- 
ues of m - from left to right: m — (brown), \m\/£ = 0.4 
(blue) and \m\/i = 0.8 (green). As \m\ changes, the Bessel 
functions get stretched and therefore for a given I the coef- 
ficients get smaller. Also plotted for comparison (black, dot- 
dashed) is the linear approximation for m = 0. Note that for 
£ 2> 5 the polarization and temperature coefficients become 
indistinguishable. 

leading order coefficients were never derived before in the 
literature. The {£, £ — 3} is for instance simply: 

6 

+ 0(I3 5 ), 

(12) 

and similarly simple expressions hold for other {£, £±n} 
correlations. Another cross-check of the above formulae 
is to confront the expansion up to second order with the 
coefficients in [5] , but there is a subtlety involved and we 
come back to this issue in Section II B. 

Figure 1 depicts some of the nonlinear coefficients (9) 
for the cases n = 1 and n — 2 and for different values of 
to. As to changes, the Bessel functions get stretched and 
therefore for a given I the coefficients get smaller. This 
plot in particular was made computing (9) for i = 100, 
but as discussed above, the coefficients are essentially 
the same for any £, as long as j3 <C 1. Figure 2 depicts 
the regions for which pi is greater than 0.2, 0.5 and 1. 
Above these regions the linear coefficients differ from the 
non-linear ones by {0.4%, 2%, 10%}, respectively. Linear 
estimates are thus correct to within 10% as long as the 
measured residual velocity lies within the /3£ < 1 region 
of this plot. This should always be the case if one pre- 
deboosts the CMB since then f3 —> (3 rcs <C 10~ 3 , unless 
the residual contribution to the CMB dipole is unexpect- 
edly large. 

One could therefore conceive of the above fitting for- 
mula as direct method to estimate the correlators which 
bypasses the need to pre-deboost the CMB. Neverthe- 
less this is a subtle point since, as shown in [7], in the 
velocity estimator there is a nearly-exact cancelation of 
the leading-order term in £ and it is thus possible that 
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Figure 2. Regions for which the first order Taylor series used 
in Section IV starts to break down, as a function of the aber- 
ration velocity. The linear coefficients differs from the full 
result by roughly {0.4%, 2%, 10%} for pi = {0.2, 0.5, 1}. If 
one pre-deboosts the CMB, one expects v TCS <C 100 km/s. 

small corrections to the value of the coefficients lead to 
moderate corrections to the estimators. One would there- 
fore need to carefully to check whether the accuracy pro- 
vided by the fits here proposed is enough for such anal- 
ysis. Moreover, as stated in Section I the usefulness of 
such method would also be dependent on whether also 
in this case three simple estimators can be written ex- 
plicitly for the three cartesian components of the veloc- 
ity, thus avoiding the need to compute and minimize a 
X 2 for all possible sky directions. If both issues can be 
circumvented, then this would provide a technique to 
measure our velocity, which could be used as a cross- 
check for the pre-deboost technique. In any case in the 
present paper we propose to use the easier route of using 
the pre-deboost technique, which is more straightforward 
and without complications, as we will show in section III. 
However, before doing that we analyze further the conse- 
quences of the approximations we have found using the 
Bessel functions. 

As a side note, the fit (9) gets worse for high values of 
P and in fact breaks down whenever P is close to 1 . This 
can be understood as one must have on physical grounds 
that P < 1 and from (2) that the integrals must approach 
zero in the limit P — > 1; the proposed lit instead does not 
go to zero in the same limit. 3 This fact is not important 
for the realistic case, since P ~ 10~ 3 , but it might be 
relevant if one tries to use our fitting function to try and 
guess an exact analytic solution to the integrals (2). In 
fact, the very high precision and breadth of applicability 
of an arguably simple fitting function (9) (which holds 



3 We found a correction to the proposed fits for /3 > 0.1: the fits 
get much better as ft — > 1 (but not exact) if one replaces ft by 
/3/(l — /3 2 ) 1 / 4 in the argument of the Bessel functions. Note that 
in practice, however, such a correction is irrelevant as we know 
that /3 ~ 10 -3 . 
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for any n and either temperature or polarization) hints 
to the existence of an analytical solution involving these 
Bessel functions, which perhaps can only be derived in 
some special limit. An analytical expression for (2) or (5), 
if it exists, could shed some light into the mathematics of 
this problem and into the limits of validity of the approx- 
imation needed to make such a derivation. Nevertheless, 
such analysis is beyond the scope of this work. 



B. Consequences of the exact coefficients 

The aberration kernel was computed up to second or- 
der in [5] and it is interesting to compare their results 
with ours. In fact, those analytic perturbative results 
motivated in part the fits (9) and (10). Note however 
that the coefficients in [5] differ slightly from the ones 
here and in [6, 7] due to the fact that theirs were ob- 
tained assuming a transformation in total intensity and 
not temperature (as it should be done) and this entices 
a small correction due to the Doppler effect. In practice, 
in [5] the integrand in (2) and (5) had an incorrect ex- 
ponent to the factor 7 (1 — (3x); to wit +2 instead of —1, 
which was first realized by [6]. This correction has two 
interesting consequences. 

First, it was shown in [5] that there should be a bias 
in the measured CMB intensity power spectrum propor- 
tional to 4/3 2 ; here we find instead that for the tempera- 
ture power spectrum (the one which is physically of more 
interest) the correction is exactly zero at C(/3 2 ). Thus 
if there is such a bias, it must be at most C(/3 4 ) 4 . An- 
other interesting result in [5] is that there should be small 
0(f3/t) cross-correlations EB and TB which are usually 
absent unless parity is violated by some exotic process. 
In other words, there should be an additional term in (1), 
to wit 
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Here we find instead, by evaluating analytically (5) for 
specific values of {£, £' , m} that this cross-correlation is 
again exactly zero, at least up to C(/3 6 ). These null re- 
sults make both effects even less relevant than previously 
thought. 



C. Correlations among non-neighboring multipoles 

The relations above allow us to estimate the amount of 
signal into the correlation functions (afj^^* a f+nrn) De ~ 
tween non-neighboring multipoles n > 1, using (9) and 



4 Note however that these results are valid only for the C{ for a 
full-sky observation, while any experiment has inevitably some 
mask, implying that pseudo Cg must be used, which leads to a 
bias which is present already at C(/3), as shown in [8]. 



n = 0, using (10). Since the amount of signal between 
these multipoles are directly proportional to the coeffi- 
cients (2), a simple (and rough) estimate of the amount 
of signal in an experiment able to measure temperature 
multipoles up to a scale £ max is given by the sum in 
quadrature of the all the coefficients up to £ max . As 
shown in [7], the amplitude of the coefficients is just 
part of the story and the smoothness of the angular 
power spectrum is also relevant, and in particular for 
Cg = CV+i there is a cancelation of the leading order 
term. Nevertheless, a comparison of the values of such 
sum in quadrature between £ and £ + n, n =/= 1, and the 
sum in quadrature of the ones between £ and £+1 should 
give a hint of the relative amount of signal between these 
non-neighboring multipoles. Since there are 2^ + 1 coef- 
ficients for each £, we evaluate the quantity 
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for n > 1 or, in the case n = , instead the quantity 
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where we approximated a sum in £ by an integral and 
where oGe m is the effective m- averaged quantity defined 

as 
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which for £ > 10 effectively docs not depend on £. As 
we will discuss in more detail in Section IV, for Planck 
/3£ max ~ 2.5 whereas for some proposed future exper- 
iments f3£ max ~ 5. The above estimate tells us that 
for /^max — 2.5 (/3£ ma x — 5) the total amount of 
signal is, in comparison with the signal in {£, £ + 1}: 
around 40% (100%), for {£, £ + 2}; around 13% (70%), 
for {£, £ + 3}- around 90% (270%), for {£, £}. These 
numbers are clearly non-negligible, contrary to what was 
claimed in [6] but in agreement with results in [24]. In 
the case of the diagonal (n = 0) correlations they seem 
to be particularly large, and in fact can be higher than 
the {£, £ + 1} signal. 

The reason why it was sometimes thought that only 
the {£, £ + 1} case was measurable stems from the Tay- 
lor expansion of the Bessel functions: for small f3£, the 
{£, £} and {£, £ + 2} cases are suppressed by a factor of 
/3, and the {£, £ + 3} by a factor (3 2 . But both current 
and future experiments will reach very small scales, way 
beyond the validity of such leading-order expansion, and 
the total signal in non-neighboring multipoles will seem- 
ingly be crucial. 

Although one could argue that the above estimates are 
not refined enough to make such claims, a thorough anal- 
ysis of the full correlations in non-neighboring multipoles 
could be quite complicated. One would have to follow a 
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similar procedure as in [7] and derive the corresponding 
estimators for the velocity for at least n — and n = 2, 
taking care to isolate all the independent correlations (for 
n = 2, for instance, there is a contribution from applying 
twice the n — 1 correlations). 

The method we propose and discuss below, based on 
pre-deboosting the CMB, circumvents all the difficulties 
stressed in this Section. 



III. THE PRE-DEBOOST TECHNIQUE 

Let n be the direction of incoming light rays in a rest 
frame S at rest with respect to the CMB and n' the ob- 
served (by our CMB instrument) direction of the same 
light ray in a reference frame S' which moves with ve- 
locity (3 relative to S. The relation between the two 
directions is the aberration effect and is given by [7, 25] 
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where again 7 = (1 — /3 2 ) -1 / 2 . Besides aberration, there 
is also the Doppler effect, which relates the observed fre- 
quency v' to the emission frequency v. 



v' = z/ 7 (l + f3 ■ n) 
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The rest-frame CMB temperature field T(n) is similarly 
aberrated and we thus observe 
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The above relation is obviously invertible, and one can 
recover the original field T(n) by applying a (second) 
Lorentz boost to T'(n') with velocity —{3. 

The idea here is to use as velocity the one given by 
the measured CMB temperature dipole. In this case, the 
dipole translates trivially into a velocity /3 dip = /3 d ; p ± 
<5/3 dip for which the current uncertainty 5f3 dlp is less than 
0.3%. We thus obtain 
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This deboosted (or "deaberrated" , as it was originally 
called in [9]) temperature field should give rise to a^ m 's 
with aberration induced correlations given no longer by 
/3 dip ~ 10~ 3 but instead by a much smaller velocity, 
of order <5/3 dip ~ 3 x 10 -6 . This means that even if 
future CMB experiments are able to probe very small 
scales such as £ ~ 10000, the product 5(3 dip £ would still 
be much smaller than unity, and one could safely ana- 
lyze these correlations using only the first order Taylor 
expansion coefficients. Moreover, in this case only the 
{£, £ + 1} correlations are relevant. Therefore, by pre- 
deboosting one can rely on the velocity estimator already 



derived in [7], and one need not worry about non-linear 
corrections nor about building a more complex estima- 
tor which would take into account the other {£, £ + n} 
correlations. 

If on the other hand one detects an aberration which is 
larger than the one expected due to £>/3 d j p , this would im- 
ply that our original assumption (to wit, that the CMB 
temperature dipole is due only to a relative velocity, or 
to some other effect which is exactly degenerate with a 
velocity in all multipoles) is incorrect, and we could be 
measuring for the first time a primordial dipole contri- 
bution. The nature of such a contribution would have 
an impact on the high-£ correlations. If some type of 
perturbations are completely degenerate with a velocity 
effect also producing the same correlations at high-£, then 
they would be of course not detected by our method. In 
general since different type perturbations may lead to dif- 
ferent predictions about the high-^ correlations, then this 
would lead to the possibility to distinguish for instance 
an adiabatic mode from an isocurvature dipolar pertur- 
bation or a dipolar lensing effect or some other exotic 
contribution such as one due to vector fields. We leave 
this important question for future work. 

The idea of deboosting the CMB was first proposed 
in [9] with the goal of simplifying analysis of general 
CMB correlations. In the same paper it was pointed out 
that applying (20) directly into the reconstructed tem- 
perature maps could lead to some propagation of errors 
(and perhaps generation of spurious correlations) aris- 
ing from an interpolation of a pixelised map. The ideal 
solution is therefore to deboost the temperature field di- 
rectly in the instrument raw time-ordered data (TOD). 
In other words, one should correct the direction and fre- 
quency of each point in the TOD by applying a — /3 d jp 
deboost with n' corresponding to the precise pointing 
of the instrument (and orbital velocity, as the rotation 
of the Earth around the Sun introduces a 10% modula- 
tion to the dipole) at the original time the photons were 
collected. The corresponding pre-deboosted a^ m 's must 
only be computed a posteriori from the deboosted TOD. 



IV. APPLICATIONS TO CURRENT AND 
FUTURE CMB EXPERIMENTS 

A. Summary of CMB Experiments 

In this Section we explore the idea in more detail and 
estimate the expected signal strength and detection pos- 
sibilities in a number of CMB experiments. As it turns 
out, experiments which cover only a few percent of the 
sky are not favored to detect the proposed correlations, 
due to extra contribution from cosmic variance (see (21) 
and (33) below). These include among others: CBI [38], 
QUAD [39], QUIET [40], BICEP [41], SPIDER [42-44] 
and POLARBEAR [45]. We will focus instead on the ex- 
periments which cover a substantial fraction of the sky, 
in particular WMAP, Planck and SPT as well as some 
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Experiment 


# v bands 


10V T (if ) 


ioV(^) 


^fwhm 


/sky 


S/N 


ACBAR '08 [26] 


1 


0.9 


- 


4.8' 


1.7% 


1.0 


WMAP (9 years) [27, 28] 


5 


14 


20 


13.2' - 52.8' 


78% 


0.7 


EBEX [29] 


3 


0.33 


0.48 


8' 


1% 


0.9 


BICEP2 (2 years) [30, 31] 


1 


3.2 


4.6 


0.6' 


2% 


2.5 


Planck (30 months) [28, 32] 


7 


1.0-8.4 


1.7-14.5 


4.7' - 32.7' 


80% 


5.9 


SPT SZ [33, 34] 


3 


5.7-30 


- 


1.0' - 1.6' 


6% 


2.0 


SPTPol (3 years) [35] 


2 


1.3- 1.5 


1.9-2.1 


1.0' - 1.6' 


1.6% 


2.5 


SPTPol Wider (6 years) 


2 


2.4-2.6 


3.3-3.7 


1.0' - 1.6' 


10% 


5.2 


ACTPol Deep (1 year) [36] 


2 


0.5 - 2.2 


0.7-3.1 


1.0' - 1.4' 


0.36% 


1.4 


ACTPol Wide (1 year) [36] 


2 


2.5 - 11 


3.5-16 


1.0' - 1.4' 


10% 


4.4 


ACTPol Wider (4 years) 


2 


2.5 - 11 


3.5-16 


1.0' - 1.4' 


40% 


8.8 


COrE (4 years) [28] 


15 


0.07-9.0 


0.12 - 15.6 


2.8' - 23.3' 


80% 


14 


EPIC 4K [37] 


9 


0.08 - 0.82 


0.11 - 1.2 


2.5' - 28' 


80% 


16 


EPIC 30K [37] 


9 


0.20 - 4.4 


0.28 - 6.2 


2.5' - 28' 


80% 


13 


Ideal Exp. (up to t = 6000) 


Any 








0' 


100% 


44 



Table I. Summary of CMB experiments. The second column relates the number of frequency channels observed; #f w hm is the 
beam size diffraction limit with full width at half maximum; or is the thermodynamic temperature sensitivity per pixel; ap 
likewise for the polarization quantities E and B; / s k y indicates the fraction of the sky covered. The quoted ranges of sensitivity 
and resolution stand for the different frequency bands below 420 GHz (which we find to be the relevant ones in the cases here 
considered). In the last column we quote the computed signal-to- noise ratio, but one should note that these values depend 
somewhat on the fiducial spectra (see text). 



of the proposed future ones (SPIDER, ACTPol, SPTPol, 
Cosmic Origins Explorer - COrE and The Experimental 
Probe of Inflationary Cosmology - EPIC). To illustrate 
the difficulties faced by surveys probing only a small piece 
of the sky we will make exceptions for some of those: 
ACBAR, EBEX and BICEP2. 

Estimates in this section refer to statistical noise alone; 
care must be taken when interpreting these due to the 
presence of foregrounds and systematic noise. 

In Tables I and II we list a summary of CMB experi- 
ments. In the former we list the range of some parameters 
on the different frequency bands; in the latter we con- 
sider only the most promising subset of experiments and 
moreover, for a more detailed comparison, consider only 
the best frequency channel for measuring the spectra. In 
both tables #f w hm is the beam size diffraction limit with 
full width at half maximum (fwhm); ut is the thermo- 
dynamic temperature sensitivity per pixel (a square the 
side of which is the fwhm extent of the beam); ap like- 
wise for the polarization quantities Q and U ; £^ v \ im is 
the multipole at which the temperature noise spectrum 
equals cosmic variance; d^ax i s ^ ne multipole for which 
the TT or EE spectrum equals the instrument (statisti- 
cal) noise (i.e., S/N = 1). For both SPTPol and ACTPol 
we assume a net 9-hour per day observation time. 5 



The quantities 9 



fwhm 



cvlim ' 



/sky and ot are re- 



5 This in principle could be made higher by making observations 
also during daytime, although this extra exposure would be done 
with less sensitivity due to heating of the telescope [46] 



lated by the expression for the noise power-spectrum (see 
e.g. [47]): 



AC, = 



/sky(2^+l) 



Ci + N t 



(21) 



where the first term stands for the cosmic variance (CVi) 
and Ni for the instrumental noise (see [47-50]): 



N, = 



9fwhm°T ex P 



£(£+1) 



Q2 
lwhm 

8 In 2 



(22) 



Cosmic variance is predominant at lower and intermedi- 
ate values of £ (to wit for t < ^lim); instrumental noise 
dominates for £ > £T„ ]lm and determines £*„. Note that 

CVllIIl IllciX 

in the case of Earth experiments the estimate (22), some- 
times referred to as the Knox formula, is inaccurate for 
I < 500 due to possible atmospheric fluctuations [36]. 
Note also that we list sensitivities per pixel in /J.K/K, 
but some of the references here listed prefer to describe 
them as either noise-equivalent temperatures (NET) in 
M-^CMBV^ (for a single detector - e.g. one bolometer) 
or as sensitivity in \xK- arcmin (for a given frequency 
band, with all detectors in that band combined). To 
convert between these quantities one has to make use of 
the following relations [37]: 



a(/j,K ■ arcmin 



' 87r/ sky [NET^cmb y/s)] 10800 



^mission ( ^ ) -^detectors 



cr((iK ■ arcmin) 
2.725if #f wmn (arcmin) 



(23) 
(24) 



Experiment 


best v (GHz) 


10Vr(<#) 


lOVp(i^) 


"fwhm 


/sky 


f T 

^cvlim 


f 

L max 


i E 

L max 


WMAP (9 years) 


94 


14 


20 


13.2' 


78% 


600 


900 


5 


Planck (30 months) 


143 


1.0 


1.7 


7.2' 


80% 


1800 


2500 


1700 


ACTPol Wider (4 years) 


150 


2.5 


3.5 


1.4' 


40% 


3600 


4600 


3300 


COrE (4 years) 


225 


0.21 


0.36 


4.7' 


80% 


3100 


3700 


3000 


EPIC 4K 


220 


0.24 


0.34 


3.8' 


80% 


3400 


4300 


3300 


EPIC 30K 


220 


0.56 


0.79 


3.8' 


80% 


3000 


3900 


3000 



Table II. Similar to Table I for selected CMB experiments which we will consider in more detail. Values quoted stand for the 
optimal frequency band with respect to a combination of temperature and angular sensitivities, gt is the thermodynamic 
temperature sensitivity per pixel; op likewise for the polarization quantities E and B\ i^viim is the approximate multipole 
at which the temperature noise spectrum equals cosmic variance; is the approximate multipole for which the TT or EE 

spectrum equals the instrument (statistical) noise. Note that or , #f w hm , /sky, i^viim an d ^max are related by equation (22). 



Note that in (23) above, / s k y is the total fraction of the 
sky covered by the corresponding CMB mission before 
any possible sky cuts (usually around the galactic plane). 
Thus, / s ky = 1 for all space based surveys and we assume 
/sky = /sky for all ground based surveys here considered. 

Also worth noticing is the relation between <jt and 
<JE;B- If both linear polarization states are given equal in- 
tegration times (as is usually the case), the total number 
of photons available for the temperature measurement 
will be twice the number available for either polarization 
measurement, and one has cfe;B = v2ct [51]. This re- 
lation, however, assumes that all detectors (being cither 
coherent receivers, like HEMT amplifiers, or incoherent 
ones, like bolometers) in the experiment are sensitive to 
polarization. If that is not the case, then one has instead 
the relation 



r, -^total 

OE-B = 4/2— CT T , 

V iVpolar 



(25) 



where iVtotai stands for the total number of detectors and 
■^Vpoi < -^totai is the number of polarization detectors. 

Finally, the measurements in the different frequency 
bands are formally independent, so one can combine dif- 
ferent channels to get lower instrumental noises (cosmic 
variance is completely correlated in all frequency bands 
and therefore cannot be mitigated) and thus probe higher 
multipoles. For Planck, as an example, combining 2 
channels (143 and 217 GHz) allows us to increase ^ ax by 
around 20%. However, more important than an increase 
in £L av , the lower instrumental noise after combining all 

ill tlx ' *— » 

channels means that the correlations here addressed can 
be detected with higher S/N ratios. To wit, combining 
the different bands one gets 



N, 



comb 



#bands 



-1/2 



(26) 



where Nj; is the instrument noise spectra for the i-th 
frequency band. 

One must nevertheless be aware that different chan- 
nels are subject to largely different foregrounds and sys- 
tematics, so that actually the usable / s k y differs for each 



frequency range. In what follows, for simplicity, we will 
assume a fixed / s k y in all frequency bands for each ex- 
periment. 



B. Detectability of our Proper Motion 

In this subsection we seek to answer the following ques- 
tion: for a given CMB experiment, what is the smallest 
value of /3 rcs it can in principle detect (say, S/N = 1 or 
3)? 

Using the induced off-diagonal ai m correlations, an es- 
timator for our peculiar velocity was built in [7]. We 
summarize below the procedure. We first define the fol- 
lowing basic quantities: 



XY 
im 



n X *n Y 



(27) 



in the appropriate frame, where X, Y stands for ei- 
ther temperature (T) or one of the polarization channels 
(E,B). The useful quantities are the real part 6 of the 
above 



rXY _ ^frpXY 
Jim — 2 \ r lm 



(28) 



Given a peculiar velocity /3 r , 



Note that f?J = // 
we can predict the average value of the f^'s, which are 
given by [7]: 



(fL Y ) = cr Y 



t+imyi 



c, 



XY 



„+X r ,XY 
La 



(29) 



6 The imaginary part has zero average. Note however that in the 
presence of some primordial parity violating effect the imaginary 
parts of FfB and Fj^ could be be nonzero, already in the CMB 
rest frame. Such observables would presumably be small and 
the correction due to a boost would have a further suppression 
due to f3, so this would most likely be a negligible correction at 
large scales. However such correction due to a boost may be 
relevant at very small scales when fit ~ 1. Investigating in detail 
such possibility can be interesting in a specific parity violating 
scenario, which goes beyond the scope of our present work. 
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Figure 3. Precision in the measurement of 8/3 achieved by different experiments using the temperature T, the E mode of 



polarization the TE cross power-spectra and the combination of all these. Here we still assume /3 rcs = 
pre-deboosting). Top to bottom: ACTPol Wide (1 year); Planck (30 months); ACTPol Wider (4 years) 
years); EPIC (4K) and finally, represented by the dashed line, an ideal experiment (no noise). 



1.231 x IO"" 3 (i.e., no 
; EPIC (30K); Core (4 



The cf* coefficients are different in the case of temper- 
ature and polarization and are given by (8); to wit: 



m 3 

c em 



fires & s G ? 



(30) 



where s — for temperature and s = 2 for polarization, 
and where we made use of (7). Note that in the limit 
of fiat spectra (Cf Y = Cfi$), one has (ffj) = 0. One 
can show, using (9), that this remains true up to C(/3 3 ); 
we come back to this issue in the end of this Section. 
Therefore (jf Y ) is generically higher the more wiggled 
the power spectrum is. 

Since (f^Y) is proportional to /3 res , it is useful to define 
a related quantity 



XY\ 



A'cs (fe m ) 



(31) 



From these predictions, we can compute an estimator 

rXY, obs ; ixy \ \ I r 
Jem \Jtm / I I \ " \Ji 



E 




cXY\2 



£e£e+i 



(32) 



where (see [7, 49, 52]) 7 




Ct + Nf) 



(33) 



and the sums are in principle carried out with — £ < m < 
£, but can also be simplified to a sum on < m < £ since 
fhn = fe^-m- ^ n approximate value for the variance of 
/3 rcs can also be written as [7] 



5(3 



2 ) 



EE E 



XY\2 
Im I 



L X 



(Ji 



N 

s 

(34) 



This leads to an important conclusion: the estimate of 
5/3 does not depend on the value of f3 rcs itself. Note that 
this is physically motivated by the simple fact that the 
error in /3 is mainly due to cosmic variance which, to 
leading order, is independent on j3. Therefore, one can 
determine the achievable precision of a given experiment 
without knowledge of the residual velocity f3 rcs . Likewise, 



7 Note that in [7] there was a typo and the square root was missing 
in (33). 
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Figure 4. Similar as Figure 3 but for the precision in the 
measurement of 8v achieved by different experiments using 
a combination of TT, TE and EE channels. Here we need 
not assume any value for /3 res . Top to bottom in both plots: 
ACTPol Wide (1 year); Planck (30 months); ACTPol Wider 
(4 years); EPIC (30K); Core (4 years); EPIC (4K) and finally, 
represented by the dashed line, an ideal experiment (no noise). 



one can also determine what is the minimal value /3 res 
that a given experiment could detect (again, say, with 
S/N = 1 or 3). 

Finally, it was also shown in [7] that the estimate on the 
direction of the velocity is directly related to the one of 
the magnitude. To wit if one chooses the fiducial velocity 
as the z-axis, for small angles the error on the absolute 
value of the velocity and the error in the magnitude of 
the direction are related by 



fires 



(35) 



Note that although as stated above the achievable pre- 
cision in 6/3 does not depend on the fiducial value of 
the velocity (and thus does not depend whether one pre- 
deboosts the CMB or not), the same is not true for 69, 
so clearly the direction of the residual velocity will be 
measurable only if /3 rcs is significantly larger than 5/3 for 
a given experiment. 

Figure 3 depicts the expected precision of some of the 
CMB experiments listed in Table I as a function of the 
highest multipole £ taken into account. Note that these 
predictions depend somewhat on the assumed fiducial 
power spectra, so quoted values should not be taken at 
face value. Here we made use of CAMB [53] and took 
its default cosmological parameters as our fiducial ones 
(see Appendix B for more details) . Experiments not rep- 
resented in the figure cannot produce a detection since 
their maximum S/N is below 1 by a good margin (except 
for WMAP, for which we forecast a S/N ~ 0.7, which is 
almost borderline). Note that even the best experiments 
(Core and EPIC 4K) fall short of the (arguable) target of 
~ 1% fractional precision on the velocity. Nonetheless, 
their precision would be able to detect a non-standard 
signal which amounts to at least 8% in the CMB dipole. 
It will take an even better experiment, with sensitivities 



Figure 5. Maximum theoretical accuracy a CMB experiment 
can achieve in measuring the correlations in the different chan- 
nels. Note that unlike all proposed CMB experiments here 
analyzed, correlations in the BB channel are relevant here, 
and in fact its S/N surpasses the ET + TE case. 



and beam sizes which are respectively high and small 
enough to probe temperature and polarization multipoles 
up to I ~ 5000 to achieve the level of 2%. Note however 
that going to smaller and smaller scales the signal gets 
more and more contaminated by local sources, and if they 
are not taken properly into account this could prevent us- 
ing this part of the spectrum, because it could introduce 
an additional preferred direction in the data. 

In Figure 4 we show a similar signal to the one of the 
bottom right plot of Figure 3, but in terms of absolute 
precision in the velocity, which does not depend on the 
value of /3 res (and therefore, on whether one deboosts the 
CMB or not). The bottomlinc is that Planck could mea- 
sure our peculiar velocity v with a precision of 55km/s, 
whereas Core and EPIC 4K/30K could do the same with 
only 20 — 25km/s of error. Finally, an ideal experiment 
probing temperature and polarization multipoles up to 
£ ~ 5000 could achieve Sv ~ 8km/s. 

We have so far ignored the B-mode polarization chan- 
nel. This is due to the fact that its spectra is supposed to 
be highly suppressed with respect to the _E-mode one. In 
fact, assuming no tensor perturbation modes (and thus 
that the BB spectra is due to gravitational lensing alone) 
we find that none of the proposed experiments here con- 
sidered can detect the correlations with S/N ratios larger 
than unity. To wit, we find that Core would have its 
S/N ~ 0.6 while for Epic 4K we would have S/N ~ 0.8. 
An ideal experiment could instead have for its BB corre- 
lations a S/N > 1 if it went until £ > 2000, and S/N ~ 10 
if it went until £ = 4000. 

Figure 5 depicts the maximum theoretical accuracy a 
CMB experiment can achieve in measuring the correla- 
tions in the different channels. Note that unlike all pro- 
posed CMB experiments here analyzed, correlations in 
the BB channel are relevant here, and in fact its S/N 
surpasses the ET + TE case. One of the reasons for this 
is the fact that for very high scales (£ > 5000) the BB 
spectra is not too much smaller than the EE and TE 
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Figure 6. Comparison between (fj^) for m = evaluated to 
first and up to third order (upper two curves, superimposed) 
and assuming a very high /3 rcs = 5 x 10~ 4 . The difference 
between both curves is depicted below (light blue curve). We 
smoothed out this last curve (as it oscillates heavily around 
zero) for better clarity. Note that even for such a high residual 
velocity the error is only ~ 1% and that for \m\ > the 
relative error is even smaller. 

ones (see Figure 7). 

The results in Section II allow us to estimate the ac- 
curacy of the first order expansion we use here, as p res i 
approaches unity. A direct estimate would be obtained 
by comparing the non-linear fitting functions (8) with 
a linear expansion. However, a more careful compari- 
son has to be done at the level of the (fe m ), because of 
the leading-order cancellations, which could be disrupted 
by small changes in the coefficients. We therefore eval- 
uated, using (9), (fem) up to third order in j3 (second 
order contributions are identically zero for the {£, 1+1} 
correlations) and compared the result with the first order 
(fem) ■ Amazingly, the leading order cancelation is basi- 
cally undisturbed, and the correction to (ft m ) is found 
to be small even for high values of /3 res : it is ~ 1% for 
/3 ros = 5 x 10~ 4 . Figure 6 depicts this property of (fe m )- 



V. CONCLUSIONS 

In this work we proposed a method to overcome the 
challenges posed to the measurement of our peculiar ve- 
locity (3 through its aberration and Doppler effect in the 
CMB at very small scales, due to the breakdown of the 
perturbative calculations when £ ^> 1//3. This technique 
consists of pre-deboosting the CMB (at the level of the 
TOD) and subsequently evaluate the estimators origi- 
nally constructed in [7], which are based on a linear or- 
der Taylor expansion in /?. The expected small residual 
velocity after such deboost justifies the use of only the 
linear order terms, and greatly simplifies the analysis. In 
particular, the pre-deboost method validates analysis of 
the aberration effect constrained to the correlations be- 
tween {£, i 4- 1} only. 

Making use of the estimator of the velocity derived 



in [7] for such correlations, we investigated the precision 
with which many of the current and future CMB exper- 
iments might be able to measure our velocity v. We find 
that Planck will put an error bar Sv ~ 60 km/s, simi- 
lar to what ACTPol could achieve in only 2 years with 
a survey covering 40% of the sky (4 years of observation 
would result in Sv ~ 40 km/s). SPTPol, could also be 
competitive if it carried out a similar wide survey. Pro- 
posed future space experiments such as EPIC and Core 
could put a limit Sv < 30 km/s. Even more precise ex- 
periments, able to measure the CMB signals (and get 
systematics under control) all the way to I ~ 5000 could 
in principle achieve Sv < 10 km/s. 

Since one expects from the CMB dipole that v ~ 
370 km/s, with current (near future) experiments one 
would be able to detect a residual term of primordial 
origin which contributes to a fraction larger that 20% 
(10%) of the dipole. This is of course true unless the pri- 
mordial effect itself also induces a dipolar kernel leading 
to correlations, with exactly the same coefficients as the 
peculiar velocity we analyze here, on the I > 1 multi- 
poles. Apart from velocity, other possible components to 
the dipole and/or a dipolar kernel include intrinsic adi- 
abatic, isocurvature and lensing effects, as well as more 
exotic effects, some of which we now briefly discuss. 

An example would be given by some fundamental vec- 
tor field or perhaps magnetic fields coherent over the 
horizon scale, which could single out a preferred direc- 
tion and therefore induce a dipole. Some inhomoge- 
neous and anisotropic cosmological models are also ex- 
pected to contribute to the dipole. For instance inho- 
mogeneous gigaparsec-scale void models, embodied by 
the Lemaitre-Tolman-Bondi metric, are sometimes con- 
structed to serve as a candidate for dark energy, and 
have recently drawn considerable attention from the lit- 
erature (see, for instance, [54]). In these spherical sym- 
metric models, the observer is sometimes assumed to be 
in the center for simplicity, but any off-center displace- 
ment induces a dipole contribution. This in fact happens 
to be the tightest constrain in this off-center distance, if 
one assumes the CMB dipole to be entirely due to this 
displacement and that any additional peculiar velocity is 
actually zero [55]. Metrics which exhibit vorticity (e.g., 
Bianchi V and VII) also can contain a preferred direc- 
tion and sense, and so also introduce a dipole [56]. On 
the other hand, homogenous anisotropic metrics with- 
out vorticity, such as Kantowsky-Sachs, Bianchi I and 
III, exhibit a preferred direction but no preferred sense 
(i.e., a preferred axis, but not an arrow), so they con- 
tribute to the CMB angular power spectrum at most at 
the quadrupole level [57, 58]. Finally, it is also possible 
that more generically large-scale vector perturbations in 
FLRW with unusually large amplitude may induce such 
dipoles, as well as a primordial power spectrum which 
explicitly depends on the direction of the Fourier mode 
k instead of just k = |k| [59]. 

It is an interesting question, which deserves future in- 
vestigation, to check what kind of primordial perturba- 
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tions or other effects such as the ones discussed above 
could induce similar correlations between different £'s and 
if the produced correlations are proportional to the pro- 
duced Doppler effect. Since it is natural to expect that 
in some cases there is no such proportionality, aberra- 
tion can provide in principle an observational handle to 
distinguish between these different contributions. 

While the prc-deboost is probably the simplest way to 
proceeed, we have also proposed very precise fitting func- 
tions for the complicated integrals (2), valid also when 
I >• 1/(3. Although another method was recently pro- 
posed in [24] to compute this kernel elements, it relies 
on recursion relations between some of its elements and 
is not as simple to implement as just using the simple 
Bcssel functions (9) and (10). 

Apart from introducing a correlation between differ- 
ent atmS, it was showed in [5] that our peculiar velocity 
would also cause a bias in the measured CMB intensity 
power spectrum proportional to Af3 2 . Here we showed in- 
stead that for the temperature power spectrum (the one 
which is physically of more interest) the bias is exactly 
zero at C(/3 2 ). We also showed that the small cross- 
correlation between E and B modes computed in [5] is 
instead also zero - at least up to C(/3 6 ) - when one cor- 
rects the aberration kernel to transform as temperature 
instead of as intensity. 

With some recent claims of the existence of possi- 
ble unexpected bulk flows [17-20] it is a very interest- 
ing prospect that we will soon have a complementary 
measurement of our peculiar velocity with respect to the 
CMB. A dipolar effect has already been measured out- 
side the CMB: in supernovae [60], in the cosmic infrared 
background [61] and also marginally in X-rays [62]; it 
could also be measured in the future using the cosmic 
parallax technique [55, 63, 64], and in the optical, ra- 
dio and gamma bands. It would therefore be interesting 
to explore the interconnections and compare the future 
precision that these different approaches can achieve. 
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Figure 7. Computed angular power spectra with the param- 
eters here listed and inclusion of lensing effects due to the 
linear part of the matter power spectrum. Top to bottom: 
TT, TE (absolute value), EE and BB modes. Note that for 
very high £ the BB spectrum approaches the EE and TE 
ones 



Appendix A: Details in the Computation of the 
Non-linear Coefficients 

The integrand in (2) is an oscillating one, and care 
must be taken in the computation of the integral. The 
frequency of oscillations increase with I (more precisely 
with I — |m|), which means that for smaller scales the 
numerical challenges increase. In fact, for £ > 8 some 
standard numerical integrating procedures start to give 
incorrect results. In order to circumvent this one must 
resort to arbitrary-precision arithmetic and make use of 
higher-than-double precision (for an alternative method 
see [24]). With arbitrary-precision arithmetic we were 
able to compute some of the coefficients (2) but as £ in- 
creases so must the numerical precision, and the code 
both gets slower and memory used gets higher exponen- 
tially with precision. The smallest scale integral we com- 
puted, £ = 700, required a working precision of 300 digits, 
took around half an hour to run in a desktop computer 
for a single value of m and required roughly 5 GB of 
RAM using Mathematica. In contrast, one £ = 200 inte- 
gral can be safely computed with 120 digits and take less 
than 5 minutes in the same machine. 
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Appendix B: Details in the Computation of the 
Fiducial Spectra 

The angular power spectra used here were obtained 
using the CMB code CAMB [53], January 2011 release. 
As mentioned in the text, the predicted S/N ratios 
here quoted depend somewhat on these fiducial spectra. 
These in turn depend not only on the fiducial cosmologi- 
cal parameters used, but also on the numerical accuracy 
of the code on small scales. For the experiments here 
considered the smallest scales are roughly £ ~ 4000 but 
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we were interested also in investigating scales all the way 
to £ ~ 7000 for the case of an ideal experiment. To 
get reliable results in these scales we set in CAMB the 
Sample Boost and both Accuracy Boost parameters to 5 
(setting some of these too high seemed to produce nu- 
merical artifacts for i > 6000, such as negative EE spec- 
trum at I ~ 7500), the Accurate EE / BB flags to true 
and used £ max = 11000, fcry max = 44000 for both scalar 
and tensor modes. The cosmological fiducial parameters 



were the default ones, such as: h — 0.7, Qbh 2 = 0.0226, 
il c h 2 = 0.114, n v = rt k = 0, n = 0.96, z leion = 11 et 
cetera. We also included the corrections due to gravita- 
tional lensing due only to the linear part of the matter 
power spectrum. It turns out that a lensed angular power 
spectrum leads to overall smaller S/N in the correlations 
here studied, in particular in the polarization channels. 
The decrease in precision after summing all channels is 
roughly 20-30% depending on the experiment. Figure 7 
depicts the generated spectra. 
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